Implementation of the refined Deutsch-Jozsa algorithm on a 3-bit NMR quantum 

computer 
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We implemented the refined Deutsch-Jozsa algorithm on a 3-bit nuclear magnetic resonance quan- 
tum computer, which is the meaningful test of quantum parallelism because qubits are entangled. 
All of the balanced and constant functions were realized exactly. The results agree well with theo- 
retical predictions and clearly distinguish the balanced functions from constant functions. Efficient 
refocusing schemes were proposed for the soft 2-pulse and J-coupling and it is proved that the 
thermal equilibrium state gives the same results as the pure state for this algorithm. 

PACS numbers : 03.67.Lx, 03.67.-a, 76.90.+d 



A quantum computer which was just a theoretical con- 
cept has been realized recently by nuclear magnetic reso- 
nance (NMR). Several methods have been proposed such 
as ion trap ||§, quantum dot |,|], cavity QED ||,§, 
and Si-based nuclear spins p| to realize quantum com- 
puters but NMR || has given the most successful re- 
sults. Several quantum algorithms have been imple- 
mented by NMR quantum computer s |9H 13| among which 
the Deutsch-Jozsa (D-J) algorithm pij has been studied 
most because it is the simplest quantum algorithm that 
shows the power of a quantum computer over a classical 
one. Most of quantum algorithms, including the D-J al- 
gorithm, have been implemented only for functions of one 
and two bits. The successful implementation of a quan- 
tum algorithm depends heavily on the number of basic 
operations which increases with the number of qubits due 
to finite coherence time. Moreover, more than 2-bit oper- 
ations require more than two-body interactions which do 
not exist in nature. It is possible to avoid such interac- 
tions, though not easy, but it increases again the number 
of total basic gates and coherence may break down dur- 
ing the computation. There have been few works that 
have performed real three-bit operations llq] so far. 

The D-J algorithm determines whether an n-bit binary 
function, 



/:{0,l} r 



{0,1}, 



(1) 



is a constant function which always gives the same out- 
put, or a balanced function which gives for half of inputs 
and 1 for the remaining half. The D-J algorithm gives 
answer by only one evaluation of the function while a 
classical algorithm requires (2 n_1 + 1) evaluations in the 
worst case. The function is realized in quantum compu- 
tation by unitary operation, 

U\x)\y) = \x)\y®f(x)), (2) 

where x is an n-bit argument of the function and y is 



one-bit. If \y) is in the superposed state, (|0) 
then the result of the operation, 



U\x){ 



l >-" 1 >)-(-l)/«|x>( |0) 



ID 



V2 



V2 



|1»/V2, 



(3) 



carries information about the function encoded in the 
overall phase. If |x) is also prepared in the superposition 
of all its possible states, (|0) + |1) + • • • |2" - l))/%/2", 
by applying an n-bit Hadamard operator H to \x) = |0), 
the relative phases of the 2™ states change depending on 
/. If / is a constant function, then the relative phases 
are all same and additional application of H restores |x) 
to |0). If / is a balanced function, |x) cannot be restored 
to |0) by this operation. It is obvious that \y), being in 
the superposed state, (|0) — |l))\/2, plays a central role 
in the algorithm but it is redundant in the sense that its 
state does not change. 

This redundancy is removed in the refined D-J algo- 
rithm [JL6[ where the following unitary operator is used. 



Uf\x) = (-l)tW\x) 



(4) 



It has been shown that Uf is always reduced to a direct 
product of single-bit operators for n < 2. In this case, n 
classical computers can do the same job by simultaneous 
evaluations because qubits are never entangled. There- 
fore, meaningful tests of the D-J algorithm can occur if 
and only if n > 2. Recently, a realization of the D-J 
algorithm for n = 4 has been reported p7[ , but in that 
work, only one balanced function was evaluated and the 
corresponding Uf is reducible to a direct product of four 
single-bit operators. In this study, we investigated the re- 
fined D-J algorithm with 3-bit arguments to find out the 
pulse sequences of Uf's, and implemented the algorithm 
on an NMR quantum computer for all the functions. 

There are 8C4 = 70 balanced and two constant func- 
tions among all 3-bit binary functions. We index the 
functions with their outputs, /(0) • ■ • /(7), expressed as 



hexadecimal numbers. For example, /ie denotes the func- 
tion of which the outputs are given by /(0) • ■ • / (7) = 
00011110. Note that U fs = -Uf FF _ x , where x is a hex- 
adecimal number equal to or less than FF. The difference 
of overall phase cannot be distinguished in the experi- 
mental implementations. Therefore, there are 35 distinct 
unitary operators corresponding to the balanced func- 
tions, and one operator corresponding to the constant 
functions. Since the unitary operator corresponding to 
the constant functions, Uf 00 , is just the unity matrix, 
there are 35 non-trivial and distinct Ufs to be imple- 
mented. 

The NMR Hamiltonian of the weakly interacting three 
spin system is given by 



H 



3 



AujiL. 



3 

i<j 



(5) 



in the rotating frame, where h z is the ^-component of 
the angular momentum operator of spin i. The first term 
represents the precession of spin i about z-axis due to the 
chemical shift, Auji, and the second term the spin-spin 
interaction between spin i and j with coupling constant 
Jij. This Hamiltonian provides six unitary operators, 
hz{0) = exp[-i9I iz ] and Jy(0) = exp[-i92I tz I jz }. In 
combination with h z (9), two other operators h x {9) and 
Ii y {6) produced by rf pulses can perform any single-bit 
operations. The coupling operator Jij{0) can be used to 
make a controlled-NOT operation. The combination of 
single-bit operations and controlled-NOT operations can 
generate any unitary operations Jig ]. 

Table 1 shows the sequences of the realizable operators 
for all the 35 non-trivial distinct CZ/'s. In the table, the 
notations h , h and h were replaced by /, S and i?, re- 
spectively for convenience. Some of Uf's are irreducible 
and require three-body interaction. The sequences of re- 
alizable operators in the table were obtained by follow- 
ing a general implementation procedure using generator 
expansion |19[ . This method includes the coupling order 
reduction technique which replaces an n-body interaction 
operator for n > 2 by two-body ones. It is noticed that 
all Uf's consist of the operators of the single-spin rota- 
tions about z-axis and spin-spin interactions only. From 
now on, we call pulses corresponding to these operators 
the soft z-pulse and J-coupling, respectively. 

The balanced functions are classified into four types 
depending on the number of Jij (9) 's included in their op- 
eration sequences. It is easy to see that no qubits are en- 
tangled in type-I functions and therefore, obviously they 
are not the cases of meaningful tests. In type-II func- 
tions, only two qubits out of three are entangled. So, 
type-II functions can be said to be the stepping stones to 
meaningful tests. In type-Ill and IV functions, all three 
qubits are entangled and the functions of these types can 
be tested only by a three-bit quantum computer. There- 
fore, the realization of type-Ill and IV functions demon- 



strates true quantum parallelism. Each sequence in Ta- 
ble I is not unique for a given function but we believe that 
they are optimal ones for implementation of the refined 
D-J algorithm. 

The whole operation sequence for implementation of 
the refined D-J algorithm is given by H-Uf-H-D to be 
read from left to right. The first and second iJ's were 
realized by hard 7r/2 and — n/2 pulses about j/-axis, re- 
spectively. Since the read-out operation D can be real- 
ized by a hard 7r/2 pulse about y-axis, the second H and 
D cancel each other to make the sequence H-Uf. 

The superposed input state is generated by the 
Hadamard operation on the pure state |0). Therefore, it 
is usually necessary to convert the thermally equilibrated 
spin state into the effective pure state. In the case of the 
refined D-J algorithm, however, the thermal equilibrium 
state gives the same results with the pure state. The de- 
viation density matrix of the thermal equilibrium state, 
Pth, is approximated by 



Pth = hz + hz + hz 



(0) 



for the Hamiltonian of Eq. g, and the density matrix of 
|0), /O p , is given by 



P P = hz + hz + hz 

+ 2hzhz + %hzhz + %hzhz 
= pth + Ap. 



Ihzhzhz (7) 



The hard 7r/2 pulse for H transforms terms of pth into 
single-quantum coherence and terms of Ap into multiple- 
quantum coherence pQ). Since the sequences for 17/ 's 
consist of only the soft z-pulse(s) and J-coupling(s) which 
are dependent only on the z-components of spin angular 
momentums, Ufs do not change the order of quantum 
coherence. As single-quantum coherence is only observ- 
able, pth and p v give the same results for this case. In 
general, the thermal equilibrium state gives the same re- 
sults with the pure state if the operation sequence after 
the first Hadamard operator does not change the order 
of quantum coherence. 

The soft z-pulse and J-coupling were implemented 
by the time evolution under the Hamiltonian of Eq. o 
with refocusing 7r-pulses applied at suitable times during 
the evolution period. Since the refocusing 7r-pulse has 
the effect of time reversal, it can be used to make one 
term in the Hamiltonian evolve while the other terms 
"freeze" pl|-p3[. We optimized this refocusing scheme 
as illustrated in Fig. U] which shows the soft z-pulse on 
spin 1 and J-coupling between spin 1 and 2 as examples. 
The evolution time, T, is 9/Auji for the soft z-pulse and 
9/(wJij) for the J-coupling. Previous schemes divide the 
evolution period into eight periods and require six pulses, 
or suffer from TSETSE effect HHt because soft pulses 
exciting more than one but not all spins were used. Since 
the difficulty of experiment increases exponentially with 



increasing number of pulses, especially soft pulses, our 
scheme greatly enhances the possibility of successful im- 
plementation. Axes of successive 7r-pulses were chosen in 
the way to cancel imperfections of pulses. For example, 
four 7r-pulses in Fig. |l](a) were applied along x, —X, —x, 
and x-axes, respectively. 

In our experiment, 13 C nuclear spins of 99% carbon- 
13 labeled alanine (CH 3 CH(NH 2 )C0 2 H) in D 2 solvent 
were used as qubits. NMR signals were measured by us- 
ing a Bruker DRX300 spectrometer. The chemical shifts 
of three different carbon spins are about 5670, —3780, 
and —6380 Hz, and coupling constants J12, J23, and 
J13 are 54.06, 34.86, and 1.03 Hz, respectively. Protons 
were decoupled during the whole experiments. Gaussian 
shaped soft 7r-pulses were 2 ms in length and hard pulses 
were about a few microsecond. The length of the total 
pulse sequence was about 600 ms in the worst case. 

We implemented all the 35 balanced and one constant 
functions exactly. Fig g shows the results for the four 
functions belonging to different types shown in Table |. 
The lines of the spectra for the remaining functions also 
indicate as clearly as ones in the figure whether they are 
positive or negative. The balanced functions are distin- 
guished from the constant function because some of the 
lines are negative. The peaks of spin 1 and 3 show up as 
doublets in Fig 0(a), (b) and (c) while that of spin 2 is 
quartet because J13 is very small compared to J12 and 
J23. Fig @(d) shows, however, that the peaks of spin 1 
and 3 are in fact quartet also. They look dispersive dou- 
blets because the neighboring lines split a little by J13 
have different signs. These results agree well with the 
theoretical predictions obtained from 



operation sequences best for our implementation were 
found using generator expansion. Experimental pulse se- 
quences were made as simple as possible by using the 
thermal equilibrium state and the new refocusing scheme. 



Tr(e mt / h p e- mt ' h I+), 



(8) 



where p is the density matrix transformed by the opera- 
tion sequence H-Uf from p t h and i+ = I x + il y - 

In the implementation of the soft z-pulses and J- 
couplings shown in Fig. |lj, the end of pulse sequence 
(t = T) can be clearly defined for the J-coupling but not 
for the soft z-pulse, because the last pulse of soft z-pulse 
is a soft pulse which is much longer than a hard pulse. 
Therefore, whole pulse sequence was arranged to finish 
with the J-coupling. Our refocusing scheme decreases the 
length of the total pulse sequence and therefore, reduces 
signal decay due to decoherence. Imperfection of soft 
pulses is thought to be the main source of the phase er- 
ror and the decay of signal amplitude of some lines. This 
imperfection is more serious in the J-coupling than in the 
soft z-pulse because out-of-phase multiplets are produced 
in the former while in-phase multiplets are produced in 
the latter. Therefore, it is very important to calibrate 
soft pulses exactly especially for long sequences. 

In summary, we implemented the complete refined D-J 
algorithm with 3-bit arguments which involves entangle- 
ment. All the operations were realized by the time evo- 
lution under Hamiltonian with refocusing 7r-pulscs. The 
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FIG. 1. Refocusing schemes for (a) I z (6) and (b) Ji2(9). 
Short and long bars represent soft and hard 7r-pulses, respec- 
tively. The angle 9 can be changed by adjusting the length of 
evolution time, T. 



FIG. 2. Observed signals for (a) f&$, (b) /s6, (c) f^, and 
(d) /id. The x-axis represents frequency increasing from right 
to left. 



TABLE I. Sequences of realizable operators for Uf's corresponding to 35 balanced functions. 



Type-I 


Ae 


S,(7r)7,(-f)fl,(-i)Ji3(f) 


/3A 


5 z (f)i? z (-f)j 12 (f)j 13 (f) 


/of 


J,(tt) 


faa 


&W/»(f)fl.(-5)Ji3(i) 


/53 


S z (|)77 z (f)J 12 (-f)J 13 (f) 


/33 


5* z (tt) 


fs3 


S,(7T)J,(-f)Ji,(f)Jl3(i) 


,/35 


S z (i)77 z (i)J 12 (f)J 13 (-i) 


/55 


«,(«■) 


Ac 


S,(jr)7,(f)fl,(f)Ji3(i) 


,/5C 


S z (-|)77 z (f)J 12 (f)J 13 (f) 


/3C 


/,(ir)5,(jr) 


/B6 


i? z (7T)7 z (-f)5 z (-f)J 12 (f) 


Ae 


7 z (f) J R z (-|)J 12 (f)J 23 (f) 


/66 


&(ir)fl,(ir) 


As 


fl»(T)7*(i)S'.(-5)Jl2(i) 


At 


7 z (f) J R z (f)J 12 (-f)J 23 (f) 


/SA 


-R z (7r)J z (7r) 


fes 


^(7T)/ z (-f)S z (f)Jl 2 (f) 


Ad 


7 z (f) J R z (f)J 12 (f)J 23 (-f) 


A9 


I z (iv)S z (Tt)R z (ir) 


Aa 


i? z (7T)7 z (f)S z (f)J 12 (f) 


/T4 


7*(-f )7i z (f )7i 2 (f )7 23 (| ) 


Type-II 


Type-Ill 


Type-IV 


/lE 


/ z (7r)S z (-f)i? z (-f)J 23 (f) 


Ae 


7,(7r)5,(-f)J 23 (f)Ji3(f) 


At 


5 z (7T)J 12 (f)J 23 (f)J 13 (-f) 


Ad 


/.(7T)5,(5)Ji,(-5)J 2 3(5) 


/l3 


7 2 (f)5 z (f)J 23 (-f)Ji3(|) 


Ab 


5 z (7T)J 12 (5)J 23 (-f)Jl 3 (f) 


/3B 


7 z (7r)S z (-f)# z (f)J 23 (§) 


/27 


7 z (f )S z (f )J 23 (|)Ji3(-f ) 


Ad 


&(7r) Jl 2 (-|) J 23 (-|)Jl 3 (-|) 


As 


/ z (7T)S z (f)i? z (f)J 23 (f) 


/r-2 


7 z (-f )5 z (f )J 23 (f)Ji3(f ) 


Ai 


5,(7r)J 12 (-f)J 23 (f)J 13 (f) 
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